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Abstract 

The impression gained from the literature published to date is that the spectrum of the 
stadium billiard can be adequately described, semiclassically, by the Gutzwiller periodic 
orbit trace formula together with a modified treatment of the marginally stable family of 
bouncing ball orbits. I show that this belief is erroneous. The Gutzwiller trace formula is 
not applicable for the phase space dynamics near the bouncing ball orbits. Unstable peri- 
odic orbits close to the marginally stable family in phase space cannot be treated as isolated 
stationary phase points when approximating the trace of the Green function. Semiclassical 
contributions to the trace show an h - dependent transition from hard chaos to integrable 
behavior for trajectories approaching the bouncing ball orbits. A whole region in phase space 
surrounding the marginal stable family acts, semiclassically, like a stable island with bound- 
aries being explicitly ?i-dependent. The localized bouncing ball states found in the billiard 
derive from this semiclassically stable island. The bouncing ball orbits themselves, however, 
do not contribute to individual eigenvalues in the spectrum. An EBK-like quantization of 
the regular bouncing ball eigenstates in the stadium can be derived. The stadium billiard is 
thus an ideal model for studying the influence of almost regular dynamics near marginally 
stable boundaries on quantum mechanics. This behavior is generically found at the border 
of classically stable islands in systems with a mixed phase space structure. 

PACS numbers: 05.45, 03.65.Sq 



1 Introduction 

The derivation of semiclassical periodic orbit formulas for the trace of the quantum Green 
function led to a deeper understanding of the influence of classical dynamics on quantum spectra. 
Closed periodic orbit expressions have been given by Gutzwiller ||^] and Balian and Bloch |^ 
for "hard chaos" systems and by Berry and Tabor ||3|, ^ for integrable dynamics. Integrability 
and hard chaos represent the two extremes on the scale of possible Hamiltonian dynamics. The 
term "hard chaos" introduced by Gutzwiller is, however, not well defined. It implies, that 
all periodic orbits are unstable and "sufficiently" isolated to allow for the stationary phase 
approximations in the derivation of the trace formula. Hence, for lack of a better definition, one 
might say that a system exhibits "hard chaos" if the Gutzwiller trace formula as it stands is the 
leading term in a semiclassical expansion in h. 
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Neither integrability nor hard chaos is generic in low dimensional bounded Hamiltonian 
dynamics. The variety of classical systems proposed as testing models for the Gutzwiller trace 
formula in the last decade exemplify the problem of finding ideal chaos in the sense described 
above. Typical Hamiltonian systems show a mixture of the two extremes. Chaotic regions in 
phase space (containing unstable periodic orbits only) are interspersed by stable islands, which 
themselves may have a complicated inner structure consisting of chaotic bands and invariant tori. 
The "stable regions" reflect locally almost integrable behavior in phase space and can be treated 
semiclassically by an Einstein-Brillouin-Keller (EBK) approximation ||, Semiclassical 
contributions to the trace of the Green function from periodic orbits inside the stable islands 
are of modified Berry-Tabor-type as has been discussed in 0, 10, 11]. 



Much less is known about a semiclassical treatment of the transition region from the stable 
island to the outer chaotic neighborhood. The outer boundary acts as a classically impenetrable 
wall in two degrees of freedom. A coupling between dynamically separated regions can only 



be described by complex trajectories |2|, 13|. Although the outer chaotic region contains 



(per definition) only unstable periodic orbits, it shows almost regular behavior near the stable 
components. The characteristics for this kind of regularity, also called intermittency in the chaos 
literature []T4| , is a vanishing Liapunov exponent Xp for unstable periodic orbits approaching the 
island, i.e. 

_ logAp 

P~ T ' 
Ip lp—*oo 

where Ap denotes the largest eigenvalue of the stability matrix (describing the linearized dy- 
namics in the neighborhood of the periodic orbit) and Tp is the period. 

I will show that intermittency in bounded systems leads to leading order corrections in the 
Gutzwiller trace formula. Semiclassical contributions from (unstable) periodic orbits approach- 
ing the boundary of a stable island show a 7i-dependent transition from Gutzwiller to Berry - 
Tabor like behavior. This allows one to extend the concept of EBK-quantization from a stable 
region over the marginal stable boundary into its "chaotic" neighborhood. 

We will derive these results for a specific example, the quarter stadium billiard, (see Fig. 

This might appear surprising on first sight, because this billiard has been introduced as an 
example of a dynamical system which is completely ergodic with positive Liapunov exponent 
[^. The stadium billiard has thus been regarded as an ideal model for a "chaotic" system. 
First indications for the validity of the random matrix conjecture relating the level statistics of 
individual "chaotic" systems to those of an ensemble of random Gaussian orthogonal matrices 
have been found in this billiard |16, 17, 18 [. Also the discovery of scarred wave functions 



[ |19f| was first made in this system. From then on, the classical and quantum aspects of the 
stadium billiard were studied intensively both theoretically and in microwave experiments pO| , 
^], making it a standard model in quantum chaology. The billiard dynamics possesses, however, 
regularities, which prevents it from being a hard chaos system. The so-called bouncing ball 
orbits, i.e. the continuous family of periodic orbits running back and forth in the rectangle, are 
marginally stable. Another peculiarity of the stadium is the existence of whispering gallery orbits 
accumulating at the boundary of the billiard. These orbits have a strong, non-generic influence 
on the spectral statistics which causes deviations from the GOE-predictions |^l|. There also 
exist quantum eigenstates, which are localized along the bouncing ball orbits (the bouncing ball 
states), or along the boundary (the whispering gallery states), in a much more pronounced way 
than the scars found along other periodic orbits. 

As a consequence, the stadium billiard is not an example of a hard chaos billiard, but may 
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serve as an ideal model for systems exhibiting both regular and chaotic dynamics. Though there 
is no stable island, the bouncing ball orbits act as the boundary of a torus enclosing an area 
with zero volume in phase space. This allows us to study pure boundary effects without dealing 
with the inner structure of the island. 

Semiclassical results obtained for hard chaos systems suffer considerable changes when deal- 
ing with intermittency as will be shown here. First of all, the contributions of the marginal 
stable family to the trace of the Green function need a special treatment [p2| , [2^ ] . This is true, 
however, also for unstable periodic orbits close to the bouncing ball family. They can no longer 
be treated as isolated stationary phase points when approximating the trace of the Green func- 
tion at finite h. As a consequence, the Gutzwiller periodic orbit formula is not valid for the 
dynamics in a whole phase space volume surrounding the marginal stable family. The size of 
this volume is explicitly h dependent and shrinks to zero in the semiclassical limit. 

Furthermore, standard semiclassical arguments referring to a cut-off in periodic orbit sums 
at a period given by the Heisenberg time are not valid for intermittent systems. Intermittency is 
related to a loss of internal time scales in the dynamics and infinitely long orbits thus contribute 
dominantly to the semiclassical zeta function at finite h. In particular, the regular bouncing 
ball states follow scaling laws for a cut-off in the orbit summation different from the estimate 



given by the Berry - Keating resummation technique |24]. 



The article is organized as follows: in section |2|, the contribution of the bouncing ball family to 
the trace of the Green function is discussed. The derivation follows mainly the ideas developed in 
[ p^ , although here we do not deal with diffraction effects. It will be shown that the bouncing 
ball orbits give no contribution to individual eigenvalues of the stadium billiard. In section ^, 
the Bogomolny transfer operator method [26| is introduced. A Poincare surface of section is 



defined, which allows us to study the near bouncing ball dynamics, but explicitly excludes the 
bouncing ball family. The Poincare map in the bouncing ball limit is derived, which allows us 
to obtain the leading contributions to the transfer operator for the bouncing ball spectrum. An 
approximate EBK-quantization of the bouncing ball states can be derived. An analysis of the 
trace of the transfer operator unveils the breakdown of the Gutzwiller periodic orbit formula for 
unstable periodic orbits close to the bouncing ball family. In section ^, we develop the concept 
of a "semiclassical island of stability" surrounding the marginally stable bouncing ball family in 
phase space and discuss the semiclassical limit of the bouncing ball spectrum. 



2 The bouncing ball orbits 

In this section, the bouncing ball contribution to the quantum spectral determinant 

D{E) = exp / dE'TvG{E') = Y[{E - En) (1) 

•^0 n 

is derived. The {En} denote the real quantum eigenvalues of the system. Following Ref. |22, 23|, 
the Green function for the quarter stadium billiard. Fig. is divided into a bouncing ball part 
and the rest, 

G{q,q',k)=Gbb{q,q',k)+Gr{q,q',k), (2) 



where k = \/2'mE/h. The bouncing ball contribution is restricted to g, q' inside the rectangle and 
contains in semiclassical approximation all paths from q to q' not reaching the circle boundary 
or bouncing off the vertical line in the rectangle. The classical paths correspond essentially to 
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free motion, and the bouncing ball Green function Gbb can be written as an infinite sum over 
Hankel functions (plus phases due to hard wall reflections). The trace then has the form |22| 



7 hi' 

TrGbbik) = gbbik) = -i^k + iabk J] H'^\2bkn) + lim ^N^{ek), (3) 

n=l 

r(l) 



where a is the length of the rectangle, b denotes the radius of the circle and is the Hankel 
function of the first kind. The real part of the trace diverges due to the logarithmic singularity of 
the Neumann function Nq at the origin. This kind of divergence is a well known phenomenon for 
the trace of Green functions and can be regularised e.g. by defining gbb{E) = gbb{E) — ^e[gbb{EQ)\ 
for Eq fixed 12^, |2|. Note, that the same regularisation techniques have to be applied for the 



true quantum trace g{E) = J2niE — En) ^ and are thus not an artifact of the semiclassical 
approximation. We now move directly to the integrated trace, which reads 



Igbb{k)= dk'~gbb{k') = Bbb{k')-i^Nbb{k)-i'^Y.-Hi\'^^^'^)- (4) 




Bbb{k') - inNbbik) - - Y: ^^e^*'"""*'"- (5) 

n=l 

and the asymptotic form of the Hankel function has been used in the last step. The sum in 
(I) is of the same kind as the contribution of a periodic orbit family in an integrable system 
(here the rectangle) as derived by Berry and Tabor |^, . The real part of the non-oscillating 
contribution is a function of k"^ only, 

Bbb{k^) = ^k\\ogk^-l) + ^^Pk^ + l^^. (6) 

The real parameter (3 originates from the regularisation and has the form 
lo£ k °° 

^ = _ y ]v(i)(26A;on) with /cq = V2mEoh. (7) 

TT 

n=l 

The imaginary part is the bouncing ball contribution to the mean level staircase function 

The term Bbbik"^) is of the general form as derived in [^8| for arbitrary billiards with compact 
domain, i.e. 

B(k'^) = —k'^riogk'^ - 1) + log k'^ + c, (9) 

An 

where A corresponds to the volume of the billiard and /3, 7, c are real constants. Note, that the 
logarithmic singularity for /c — > in (P) is absent in (P). The sum in (Q) is convergent on the 
whole complex fc-plane (after analytic continuation for Im{k) < 0). The function Igbb{k) is free 
of singularities and analytic everywhere apart from the lines Re{k) = rmr/b, m integer, where 
it exhibits square - root cusps. The dominant oscillation introduced through the bouncing ball 



orbits is clearly visible in the oscillatory part of the level staircase function |21, 23 1, and agrees 
with the prediction (^. The bouncing ball part does, however, not contribute to individual 
eigenvalues. The spectral determinant (|l|) can be factorised according to 

D{k) K Dbb{k)Dr{k) with Dbb{k) = e^3''^^^\ (10) 
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Figure 1: The modulus of the boundary integral spectral determinant before (a) and after (b) 
dividing out the bouncing ball part, Eqn. ([To|), here for a = 5, 6 = 1; (in addition Dhb{k) is 
plotted in (a)). The k - interval shown contains about 370 energy levels, which are the real 
zeros of the spectral determinant. 
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and the bouncing ball contribution Di)b{k) is a function which has no zeros in the whole complex 
k plane. It gives rise to oscillations in the modulus of the spectral determinant with an amplitude 

growing like exp I f I and a period ir/b, but it does not influence the individual zeros. The 
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bouncing ball contributions to the spectral determinant can thus be divided out without losing 
information about the spectrum. This can be demonstrated by considering the boundary integral 
determinant 

Db^mik)=detil-K{q,q',k)), (11) 

where K is the boundary integral kernel defined for q and q' on the boundary of the billiard, 
(for details see |^0|). The function Diii^ik) has the same zeros as the exact spectral determinant 
D{k), but may differ otherwise. Especially, a smooth part exp(i3(A;^) — i'irN{k)) is absent in the 
boundary integral determinant. This term originates from the zero length limit in the Green 
function and is thus a pure volume effect. The modulation due to the bouncing ball contribution 
is clearly present in Dhim{k), see Fig. |l]a. The large overall oscillation of the amplitude covers 
7 order of magnitudes for A; ~ 30 but is completely removed by factorising out the bouncing 
ball part Di,i){k), see fig. P?. Note, that the boundary integral determinant after factorization 
exhibits cusps at k = rmr/b, and is thus non-analytic there. 



As a somewhat paradoxical result, we obtain that the bouncing ball family is responsible for 
strong modulations in the trace as well as in the spectral determinant. It does not, however, 
contribute to individual eigenvalues, and does not explain the regular modulation in the level 
spacing itself. The bouncing ball family covers a volume of measure zero in phase space, and 
is not sufficient to form a support for an eigenfunction alone. (This is of course true for any 
invariant family e.g. in an integrable system.) As a consequence, the neighborhood of the 
bouncing ball family in phase space must be responsible for the existence of the bouncing ball 
states and for the periodic change in the level spacings leading to the modulation in the level 
staircase function. The next section is devoted to a proper semiclassical analysis of the near 
bouncing ball dynamics. 



3 The semiclassical transfer operator 

So far, it has been widely believed, that the remainder term in the trace of the Green function 
(^) as well as in the spectral determinant ( [To| ) can in leading semiclassical approximation be 
described by the Gutzwiller periodic orbit formula [||, or equivalently by the semiclassical 
spectral determinant [p^ j. This belief is confirmed by studies of the Fourier transform of the 
(smoothed) quantum spectrum, which agrees reasonably well with periodic orbit predictions 
[2^ , |23|| . The semiclassical expression for the determinant, valid for systems with unstable, 
isolated periodic orbits only, has the form 

D{E) = A{E)e-^(^^ exp f - V V -M^rSAE)/n - ^ra,./2)\ ^^^^ 

\ p r=i rJ|det(l-M-(^)p)| / 



The sum is taken here over all (single repeats of) periodic orbits of the system. The classical 
action Sp = § pdq is taken along the periodic orbit, and M denotes the reduced Monodromy 
matrix, which describes the linearized dynamics in phase space perpendicular to the periodic 
orbit on the energy manifold. The Maslov index a counts twice the full rotations of the (real) 
eigenvectors of M around the orbit (plus two times the number of hard wall reflections). The 
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prefactors in front of the periodic orbit product are due to the zero-length limit in the Green 
function in a similar way as derived in (P), and do not contribute to the spectrum. Eqn. (12) 
is formal in the sense, that it is not convergent for real energies [31] and suitable resummation 
techniques have to be applied 24, 33, 34 1. Thereby, the exponential function containing 
the periodic orbit sum is expanded and the resulting terms are regrouped by ordering them 
with respect to the total action or the total symbol length (after choosing a suitable symbolic 
dynamics). Such techniques have been shown to work successfully for hard chaos systems. 

In the following, I will show that formula (12) is not valid for periodic orbits within a phase 
space region surrounding the bouncing ball family and will give explicit bounds for this area. 



3.1 The Bogomolny Transfer Operator 

Our starting point is the Bogomolny transfer operator |26] which is a semiclassical propagator 
for a classical Poincare map. It has the form 



T{q,q',E) 



(27rin)(/-i)/2 



E 

cl.trq—rq' 



dqdq' 



^iS{q,q' ,E)/n-ivK/2 



(13) 



where g, q' are points on an appropriate Poincare surface of section in coordinate space and 
/ denotes the number of degree of freedom. The sum has to be taken over all classical paths 
from q to g' crossing the Poincare surface only once with momentum pointing in the direction 
of the normal to the surface. Again, S{q,q';E) denotes the classical action along the path for 
fixed energy E. For billiard systems, it equals k times the length of the classical trajectory. 
The integer number u counts the number of caustics in g-space (plus again twice the number of 
hard wall reflections). The semiclassical eigenvalues are given by the zeros of the determinant 
det(l — T{q,q' ,E)). Evaluating the determinant in a cumulant expansion using the Plemelj- 
Smithies formula [E6|, ISTtl 



det(l - T) 



E(-i: 

n=0 

1-TrT 



n 



1 



(TrT^ - (TrT) 



(14) 



with 



an(T) 



TrT 
TrT^ 
TrT^ 



n — 1 
TrT 
TrT^ 





n - 2 
TrT 



TrT"^ TrT" "'^ ^YtT^^'^ 







TrT 



(15) 



provides the connection to the expanded periodic orbit formula (^) using the iterates of the 
map as an expansion parameter. Periodic orbits appear as stationary phase points in the various 
traces and the amplitudes are recovered using the relation 





d^S{q,q') 




V|det(l-M)| V 


dqdq' 


q=q' 



d'S{q,q) 



dq^ 



(16) 



The stationary phase approximation demands periodic orbits to be unstable and sufficiently 
isolated. The last condition will be discussed in detail later. Note that the determinant and 
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Figure 2: The quarter stadium billiard; the trajectory shown corresponds to {m,l) = (—1,1). 



the cumulant expansion (14) are well defined only if the operator T is trace class, which means 



essentially, that the trace of T exists and is finite in any basis. (For further details see |37|, 
other expansion methods using Fredholm theory are discussed in [^]). 

The Bogomolny transfer matrix method has been shown to work satisfactorily for hard 
chaos 



mixed ^] as well as integrable |^, ^ systems, and has also been applied 
successfully to the stadium taking the billiard boundary as the Poincare surface of section . 
We seek now a Poincare map which reflects the whole dynamics but excludes the bouncing ball 



orbits in order to obtain a semiclassical expression for the remainder Dj. in (10). Choosing the 
intersection between the rectangle and the circle (the vertical dotted line in fig. ^) fulfills this 
criteria. The transfer operator corresponding to this Poincare section takes on a particularly 
simple form for both the near bouncing ball limit and the whispering gallery limit. It is thus 
preferable to other possible choices such as e.g. the circle boundary. 



3.2 The classical Poincare Map 

In the following, I will restrict attention to the stadium with 6 = 1. The spectrum for general b is 
obtained by simple scaling relations. I will consider the classical Poincare map (x, (p) — > {x', (p') 
with X S [0, 1] being the coordinate on the Poincare plane starting from the bottom line. The 
angle (j) £] — 7r/2,7r/2[ corresponds to the momentum vector pointing away from the circle 
measured in the clockwise direction, see Fig. ^. (The corresponding energy dependent area 
preserving map is obtained in the coordinates {x,px) = {x, \j2mE sin (/>).) The map can be 
written as 



(-1 

(-1)'"0 
arcsin(x cos 



a; + 2a tan ( 



(m + 
with 



1 



-11 



-) 



(17) 



m < X + 2a tan (/) < m + 1 



■ + 2(/ + l)6l-A;7r with l< 



+ 



TT ■ 



2Q 



< / + 1 



Figure 3: Members of the periodic orbit families approaching the bouncing bah orbits: a) (m, 
I) = (10, 0) and b) (m, I) = (-11, 1) 



The coordinates {x, (jj) correspond to the first return at the Poincare plane with momentum 
pointing toward the circle. The angle 6 is the angle of incidence for reflections on the circular 
section of the boundary; see Fig. |2[ 

The length of a trajectory for one iteration of the map is 

2a + cos(0+f) cos{^-e) 

L{x,4>) = ] h2/cos6'H — . (18) 

cos (p cos (p' 

The general expression for the Monodromy matrix is given in Appendix A. The integer numbers 
m € [—00,00] correspond to |m| -reflections on the bottom or top line in the rectangle, the sign 
of m equals the sign of (j). The index / G [0, 00] counts the number of free flights in the circle 
and corresponds to (/ + 1) - bounces with the circle boundary. The total number of bounces 
ritot with the billiard boundary is thus 

ntot = \m\ + (/ + 1) + 2. (19) 

The map provides automatically a symbolic coding. It should be noted that the code does 
not form a "good" symbolic dynamics in the sense of Markov partition theory. Multiple iterates 
of the map are not uniquely encoded by symbol strings . . . , (m, l)i, (m, . . ., i.e. the partition 
is not generating. In addition, there is strong pruning, i.e. some of the possible symbols strings 
are not realized by trajectories of the map. The symbols defined by the map (^) do, however, 
reflect the important contributions of the dynamics to a semiclassical description as will be seen 



in the next section. The problem of finding a "good" symbolic dynamics for the stadium |39, Q 
might be a consequence of this fundamental dilemma. 

The bouncing ball limit \m\ — > 00 can be reached only for Z = or 1, but for all x on the 
Poincare surface of section. The opposite limit / ^ 00, (the whispering gallery limit), is possible 
only for m = or 1 and a decreasing x - interval of starting points. Together, there are four 
infinite series of fixed points of the map which are present for all parameter values a. In the 
I = case, m must be even and positive. The corresponding periodic orbit family (2n, 0) , see 
Fig. has initial conditions 

X{2n,o)='^, <^(2n,o) = arctan(n/a) with n = 0, 1,... (20) 

and 



^{2n,o) = V(2™)' + (2«)' + 2, det(l - M(2„,o)) = 2L(2„,o) - 4. (21) 

Note that det(l — M) grows linearly with the length L and thus with the order parameter n 
which is in contrast to exponential growth for strictly hyperbolic dynamics. 
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The family with / = 1 approaching the bouncing ball orbits are formed by periodic orbits 
which start in a corner of the rectangle and bounce off perpendicular to the bottom line in the 
circle, see Fig. In this case, m must be odd and negative, i.e. m = —2n — 1. The angle of 
incidence for the bounce in the circle fulfills the condition 



2 sin^ 6* + i (2n cos 6* + 1) sin 61 - 1 = 0, n = 0, 1, , 
a 



(22) 



with approximate solution 



2n + 1 

The starting conditions on the Poincare surface are 
1 vr 



X 



(-2n-l,l) 



2cos0„' -(-2"-i'i) 2 
One obtains for the length of the periodic orbits 



+ 29n < 0. 



L 



(-2n-l,l) 



and for the weight 



2^ (cos 9n + n)2 + (sin 6'„ + a)2 + 2 cos On 

(2n)2 + (2a)2 + 4 - ^ + 0{n~^). 

2n^ 



det(l - M(_2„-i,i 



L 



(-2n-l,l) 

cos 9„ 



8n + 5— + 0{n- 
n 



(23) 

(24) 
(25) 

(26) 



(see also appendix Again, the determinant increases linearly with the length of the periodic 
orbits. 

An analysis of the periodic orbits (0,/) and (1,Z) approaching the whispering gallery limit 
/ ^ oo is provided in appendix |^. 

3.3 The Poincare map in the bouncing ball limit 



The map (17) can be considerably simplified both in the bouncing ball limit and in the whispering 
gallery limit. The latter is postponed to appendix 

For trajectories (m,0), one has to distinguish between four different cases: (m > 0, even), 
(m < 0, even), (m > 0, odd) and (m < 0, odd). One obtains for (m > 0, even) 



x' 


X 


1 


4a2 x(l - 


-X)2' 


1 - 2x 


^2 


1 - 


2x 




vr 2a 
2 dfji 


(1 


- 2x) 


1 - 


4 a2 


with 












dm 


= 2a tan (p 




m — X + x 





i-x{3- ^2) 
1 - 2x 



+ 0(m' 



(27) 
(28) 



The length of a trajectory is 

L{x, (j) 



{l-xf 2a2 Aa'^x'^{l-xf 2a^ 



(29) 
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with = X + 2atan(/) — m > 0. Note that x' depends at leading order on x only (and 

not independently on both x and (/>). For the transfer operator (|l^), we need the length of a 
trajectory as function of the initial and final points x and x' . One obtains 



L(x, x) 



with 



m — X + x + 2 + — — 

dm 



a 



+ 0(771" 



dl, + (2a) 2 + x' 



X 



1 + 2x 



- + 2 + 0(jn 



x,x' e [0, 1] 



dmix, x ) = m — X + 



(30) 

(31) 
(32) 



l + 2x' 

The length of a trajectory after one iteration of the map is thus given by its length in the rectangle 
plus corrections which depend on x' only (up to 0(777"^)). The mixed second derivatives of 
L{x, x') showing up in the transfer operator (^) are 



dxdx' 



+ 0(777- 



((i2^ + (2a)2)3/2(i + 2x')2 

The length spectrum for the Poincare map with 777 odd or negative is obtained by the 
following replacements in (RO) - (p3|): 



X 



(771 < 0; even) 777 — > —m 
(771 > 0; odd) 777 — > 777 + 1 
(777 < 0; odd) 777 — > —777 — 1 X 



-X X 



-X 



-x' X G [0, 1] 
-x' X £ [0, 1] 
X G [0, 1] 



x' G [0,1/3] 
x' G [0,1/3] 
x' G [0,1]. 



In a similar way, the approximate map for the (777, 1) - trajectories can be constructed. In 
the bouncing ball limit 777 — > 00 only (777 > 0, even) and (777 < 0, odd) is possible. Again, we 
discuss first the case (777 > 0, even), for which we obtain 



Ax - I 



1 + 



8a2 X 



(5^2 - 4x + 1) 



+ 0(7f7" 



Ti 2a , 
■2+^('^-') 

u>rn. 



4 2x(3 - x2) - 1 
3tE 4x - 1 



+ 0(77l" 



with 

dm. 



m — X + X = 2a tan 
The length of the trajectory as a function of the initial and final points is 



L(x, x) 

with 
and 



m — X — X + 4+ — — 

dm, 



dl, + (2a) 2 - x' 
dm{x, x) = m — X + 



a 



X 



Ax' - 1 

x' 
Ax' - 1 
4a2 



+ 0(777" 

4 + 0(777"^) 

X G [0,1]; x' G [1/3,1] 



dxdx' {dl, + (2a)2)3/2(4^' _ 1)2 • 
We recover the (777 < 0, odd) case by replacing 

(777 < 0; odd) 777 — > —777 — 1 X — > —X x G [0, 1]; x' G [1/3, 1] 
in dH) - (H). 



(34) 
(35) 



(36) 

(37) 
(38) 
(39) 
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3.4 The transfer operator T in the bouncing baU hmit 



The main advantage of a semiclassical description of quantum mechanics is to study directly 
the influence of (parts of) the classical dynamics on quantum phenomena. The importance of 
the near bouncing ball dynamics on the quantum spectrum can be understood by analyzing the 



T-operator obtained from the Poincare map in the bouncing ball limit, see section 3.3 
The T -operator in the bouncing ball limit can be written as 



T{x, x'] k) = To(x, x' , k) + 



To{-x,-x',k) if 0<x'<l/3 
Tilx,x',k) if l/3<x'<l 



(40) 



with 



To{x,x') 



Ti{x,x') 



2a\ 



ik(2+x'- „ 



-)— iItt oo 



2m 



l + 2x' 



-2a\ 



n=Q 
-)— i|7r CO 



27ri 



ix' 



E 

=0 

c 

E 



oJ/cLq (n) 



(Lo(n))3/2 (L+(n))3/2 

^ikL~ (n) ^ikL'^ (n) 



n=0 



(Lr(n))3/2 (L+(n))3/2 



(41) 



where L^^-^{n) is defined as 



L±(x,x';n) = W(2n ± x + ^^^)2 + (2a)2; 



Lf{x,x';n) = J {2n ± x + ^^^)^ + {2af 



(42) 
(43) 



Here, the lower index corresponds to / = or 1. The upper index — or + distinguishes between 
contributions originating from trajectories with m even or odd. Note that trajectories in the 
bouncing ball limit have only one caustic in g~space inside the circle before returning to the 
Poincare plane both for / = and 1. The additional phases from hard wall reflections, see (19), 
have already been incorporated. 

The operator (|40|) is the main result of this paper. It contains the dynamics in the stadium in 



a somewhat counterintuitive way. All contributions of trajectories with more than 2 reflections 
on the circle boundary are neglected. In addition, the short orbits for / = or 1 are repre- 
sented least accurately. This seems to contradict our common understanding of a semiclassical 
treatment of quantum mechanics for classically "chaotic" systems. The quantum spectrum for 
hard chaos systems is expected to be build up collectively by all unstable periodic orbits and 
the shortest periodic orbits are supposed to dominate an expansion of the spectral determinants 
(12) or (|T^). The stadium billiard posses, however, a subset of regular eigenstates, the so-called 
bouncing ball states, which show a nodal pattern very similar to the checkerboard pattern ob- 
tained for the unperturbed rectangular billiard. It is this subset of states which can be treated 
by the approximate transfer operator (|40| ) alone. This is shown in Fig. ^. The leading terms 
in the cumulant expansion (|l4|), i.e. det(l — T{k)) ~ 1 — TrT(A;), is plotted here as a function 
of k. The quantum eigenspectrum of the quarter stadium, marked by crosses on the k axis, is 
obtained from the boundary integral method. The eigenvalues having a bouncing ball nodal 
pattern are emphasized by arrows. (The bouncing ball states have been identified by inspecting 
individual wave functions). The minima of — 1 — TrT(A;) — coincide very well with the eigen- 
values corresponding to bouncing ball states found by our subjective criteria. The T-operator 
constructed from a Poincare map in the bouncing ball limit fails in other regions of the spec- 
trum. Some of the states are either completely ignored, see e.g. at k ^ 7.6, or they appear as 
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Figure 4: The absolute value of 1 — TrT for a = 5 together with the quantum eigenvalues (+) 



and the quasi EBK - solutions (A), see Eqn. (pQ) in section 3^. The bouncing ball states are 
specially marked by arrows. 



doublets, where the T - operator expects only a single state, see e.g. around k w 6.7. The latter 
case corresponds to a bouncing ball state interfering with a nearby state originating from the 
non-bouncing ball dynamics. By using higher terms in the cumulant expansion (or even the 
full determinant), non-bouncing ball behavior can partly be resolved. A quantization of the full 
spectrum cannot be expected, as important parts of the dynamics have been neglected. Note, 
that the trace of T has a cusp at Rek = mvr, m = 1, 2 . . .. 

In the next section, the influence of periodic orbits which appear as stationary phase points 
in the trace of the T-operator will be studied in more detail. 



3.5 The trace of T and periodic orbit contributions 

Let us first concentrate on the trace of the Tq - part in the transfer operator (|40| ) which contains 
contributions from trajectories with only one bounce on the circle boundary. The result for TrTi 
will be given later. Expanding the length terms L^{n) up to 0{n"'^) in the exponent and to 
leading order in the amplitude, see (^), (^), yields 

oo 

(44) 

-.1— V ^_pjfc(\/(2^iF+(2^+2)-i^- 

2 V ivr n3/2 

n=0 
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dx- 



1/3 l + 2x 



exp(^fc r, ) — exp(22«;x — ik 



with 

A=^(x) = x± 



l + 2x 

The phases in front of the trace - integrals correspond to the length of the periodic orbits in the 
family {m,l) = (2n,0), see (21). The negative region of integration comes form the second Tq 



term in (|40|). 

The dominant contribution to each integral is given by the first term containing A^ which 
derives from trajectories with an even number of reflections in the rectangle. It is stationary 
for X = 0, i.e. at the starting point of the periodic orbits (pO|). This becomes obvious after the 
change of variables, 

dx exp(zfc— ^ = 2 / dy exp(zfc^- . 45 

^1/3 1 + 22; 2n^ Jo V 1 + n 

Approximating the integral straight forward by stationary phase, i.e. shifting the limits of inte- 
gration to infinity, would lead to the standard periodic orbit amplitudes using ([T^), here for the 
periodic orbit family (m, /) = (2n, 0), see (pi]). The width of the Gaussian is, however, increasing 
with n and the finite limits of integration become important for n > a\/k/-K. The individual 
contributions to the sum ( |4^ ) thus show a /c-dependent transition, i.e. 

TrTo„(A:) ^ J_^^k(./WFTi^+2)-^^-f n ^ aJ - (46) 



«./^i^e*fc(V(2")^+(2a)^+2)-.i^ if n»a\l-. (47) 
2 V TT n.'^'^ ~ 




Contributions from short trajectories in ( |4q ) have the Gutzwiller form for isolated unstable 
periodic orbits with amplitudes |det(l — M„)|~^/^ n~^/^ (see (^T])) being independent of 
k. In the other limit n ^ a\ —, we obtain ^^o^'iP' l"^"^ ~ f within the integration boundaries 
and the stationary phase approximation is no longer applicable. All trajectories in the range 
of integration give essentially the same contribution to the trace as the periodic orbit itself. A 
whole manifold of orbits build up the semiclassical weights in the trace in the same manner as 
manifolds of periodic orbits on tori with commensurable winding numbers do in a semiclassical 
treatment of integrable systems [Q] . In our treatment as well as in the Berry-Tabor approach Q , 
the trace can be performed directly and no additional stationary phase approximation is needed 
(in contrast to the Gutzwiller formula). A comparison of (^) with the Berry-Tabor weights 
obtained for contributions in the rectangle, as e.g. in Eqn. (^), unveils the similarity. Note 
that the weights are now explicitly A;-dependent and decrease like n~^/^. The transition from 
"semiclassically integrable" to "chaotic" behavior affects mainly the amplitudes. The phase is 



in both cases essentially given by the length of the periodic orbit (21). 

The conceptually different treatment for integrable and hard chaos systems can thus be 
rediscovered when studying intermittent dynamics near marginally stable boundaries. The con- 
tributions to the trace interpolate smoothly between the two extremes, the transition region 
itself is, however, fi - (or k -) dependent. To the best of my knowledge, this has been shown 
here for the first time explicitly. The number of terms corresponding to contributions of isolated 
periodic orbits increases like \fk^ the turnover occurs for 

(48) 
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Figure 5: The modulus of TrTo,n as function of the index n for different fc-values showing a 
fc-dependent transition from "chaotic" to "integrable" behavior with increasing n; the fuh hne 
corresponds to |1 — M„)|, see (pB, expected for isolated periodic orbits, the dashed lines is the 
fc-dependent asymptotic form (p7|). 



The transition is indeed rather sharp as can be seen in Fig. ^ The modulus of TrTo,n is plotted 
here versus n for different k - values. The integral (44) has been calculated numerically using 
the full length formula (|3l|). The oscillations occurring in the transition region correspond to 
the phase change from 37r/2 to 7-7r/4 from Eqs. (|46| ) to (^). 

The contributions to the trace from the second term in the integral ( ^4|) vanish like 
for large k. They become, however, important for small k, especially for the /c-region below the 
ground state. 

In a similar way, the leading contributions to the Ti - operator in (^) (deriving now from 
the Lf part, see ( ^3|) and thus from trajectories with m < 0, odd) show a transition 



TVTi,„(fc) 



1 ^»fc(V(2n)2+(2a) 



2n2 ' 



a /A; log3 ^jfc(V(2n)^ + (2a)2+4-^)-47 

4 V vr n3/2 




if n <^ a\l — 
vr 



if n a\ 



(49) 
(50) 



The phase is again essentially given by the lengths of periodic orbits from the family (m, /) = 
((— 2n — 1), 1), see (p5|). A /c-dependent transition occurs in the amplitude as in (^), (^) at a 
critical summation index approximately given by the estimate (48). 

Our analysis suggests that a marginally stable boundary as provided here by the bouncing 
ball family is smoothly connected to the outer ergodic regions in a semiclassical treatment. We 
expect that this behavior is generic for systems with a mixed phase space structure. Stable 
islands are surrounded by a "semiclassically stable" layer and there is a smooth transition for 
semiclassical contributions from both sides of the classically disconnected regions. This effect 
has already been observed by Bohigas et al. ^] . These authors compared the true quantum 
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spectrum in the quartic oscillator with an approximate EBK quantization of a stable island in 
this system. They were able to attach quantum states to EBK results even beyond the boundary 
given by the stable islands. This indicates an effective enlargement of the stable region into its 
"chaotic" neighborhood by quantum effects. Our analysis provides a consistent semiclassical 
interpretation of this phenomenon. 

The results so far are based on the particular choice of the Poincare surface of section. 
The approach presented here is opposite to a semiclassical quantization of the system on the 
billiard boundary. The trace of the transfer operator contains then only contributions from the 
shortest periodic orbits, especially from the marginal stable family. The Bogomolny method 



works for this section as well [30|. Note, however, that corrections to the Gutzwiller formula due 
to intermittency are introduced here through non-classical trajectories, when approximating 
the various traces by stationary phase. Note also that a naive summation using the Gutzwiller 
periodic orbit weights (^^) for the stadium would give contributions of the form X^r^i ^ ' 6* " 
which lead to true poles at k = mvr. It is the n~^/^ fall off for large n that prevents the trace as 
well as the spectral determinant from diverging at these points. 

One might speculate why the corrections found here have not been observed in such a 
well-studied system like the stadium billiard. Previous studies were mainly based on Fourier 
transformation of the quantum spectrum. The Fourier transform exhibits peaks at positions 
corresponding to the length of periodic orbits. The Fourier spectrum is most sensitive to phases 
of semiclassical contributions and only the shortest periodic orbits are resolved when dealing with 
smoothed spectra or a finite energy interval. Intermittency as introduced through the bouncing 
ball family affects mainly the amplitudes and contributions of long trajectories. In addition, 
a /c-dependence of the amplitudes as in (p7|), ( |50| ) is washed out by the Fourier integration. 
Fourier transformation is thus very insensitive to the influence of marginally stable behavior. 
In the next section, I will show that the bouncing ball spectrum indeed originates from series 
of orbits approaching the bouncing ball family and that there is no natural length cut-off for 
trajectories contributing dominantly to the transfer operator. 

3.6 Discrete representation of the transfer operator and quasi EBK formulas 



For a further analysis of the operator (|40|), we proceed in a very similar way to Ref. [g^]- The 
infinite sums in (41) are convergent for Im{k) > and may be expressed as 



^ikV^{n) gifc-Lj^(O) +00 (-co ^ikL^{n)—2vimn 



S(4(-))^/^ 2(L^(0))3/2+^5oo^o '^'^ (4(n))3/2 (51) 
using Poisson summation [B^] . The integral representation (|5^) provides an analytic continuation 



of the kernel ( ^0|) for Im[k) < after rotating the axis of integration dn — > itidn |3^, 36|. 

For k real, the integrals in ( |5l| ) can be evaluated by stationary phase. The prefactors decrease 
not faster than algebraically and thus vary slowly compared to the phase. The stationary phase 
condition yields 

k — - 27rm = 0. 52 
on 



The solutions of 

_ 1 [ 2am x' ] 

2 yJk^/TT^—m^ l + 2x' 
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are real for m < k/-K only. The saddles for m < give no contribution due to the limits of 
integration in (|5l|) . Note, that the sum is not necessarily dominated by short orbits. On the 
contrary, infinitely long trajectories give the main contributions to the sums in ( |4l] ) when k 
approaches mvr from above. The singularities at k = rmr are linked to the cusps appearing in 
the bouncing ball contribution to the integrated trace (|^ when summing over all repetitions of 
the marginal stable family. 

By skipping the slowly varying first term on the LHS of (^Tj), i.e. neglecting again short orbit 
contributions, one obtains in the stationary phase approximation 

e2^*(i+Tf^)+*f , . 2 ' 

Toix,x') ^ -i- ^ e^Wfe'A^-^^e'T+a^^ sin(m7rx) (54) 

^ m=l 
e4fci(l-3^)+if oo^ ^ 

Ti(x,x') ^ i- y e2™Vfc'/^2-™%*3^sin(m7rx). 

4x' — 1 , 

m=l 

The amplitudes are now no longer A;-dependent. The sum over m disappears when writing the 
kernel (40) in the basis = -v/2 sin(7rn2;), n = 1,2, . . ., i.e. 

dx / dx'(pm{x)T{x,x')(pn{x'), 

/o Jo 

which leads to the discrete transfer matrix 



Tm,n{k) ^ e2..(aVW^+f + i) R^^^{k). (55) 

The matrix R is given as 

R^n{k) = -^t dx'e'^^'^''l^^^'-%'-Tf^'^^^^ (56) 
7-1/3 1 + 2x' 

^ dx'e2*^(^-2^"/(^^'-i))f'""'^ sin(n7rx') 



ii/3 4x' - 1 

and is independent of the billiard length a. Applying the unitary transformation 



T = U-iTU with Um,n = e-^-V^V-^-n^j^^^ ^ (57^ 
leads to the more symmetric form 

fnAk) ^ e2'^4f(^/^'7^+^/^V^^)+f+i] r^^^^j,^ (58) 
The transfer matrix ( |58[ ) is now essentially finite with an effective dimension 

dime// = [fc/vr] , (59) 



where the brackets [ ] denote the integer part of fc/vr. Formula (|59| ) corresponds exactly to 
the estimate for the dimension of the T-operator given by Bogomolny for the Poincare 
surface chosen here. Neither the dimension nor the matrix R in (^) depend on the billiard 
parameter a which enters only through the phase in (55), (58). The T-operator derived from 
the Poincare map in the bouncing ball limit is expected to reproduce best the bouncing ball 
states. The determinant can be approximated by its leading term in the cumulant expansion 
(14), i.e. det(l — T) « 1 — TrT. The phases in front of the matrix R in (^) yield a quantization 
condition for the bouncing ball spectrum. The matrix R acts as a filter determining the k — 
intervals, which allow for bouncing ball states in principle. A closer analysis shows that the 
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diagonal elements Rmm are dominant and approximately real in the region k £ [mir, {m + 1)tt]. 
This leads to an EBK-like quantization condition 




M^ + - + \ = M + N for M = 1,2,... (60) 
vr 4 



N = 1,2, ...,[aV2M + 1 + 5/4]. 

M and N act as approximate quantum numbers corresponding to bouncing ball eigenfunctions 
with (M — 1) nodal lines perpendicular and {N — 1) nodes parallel to the Poincare surface. 
(An additional M on the RHS of (^) is introduced for convenience). The cut-off in the - 
quantum number originates from the fe-window given by the R matrix. The states with fixed 
quantum number M are restricted to the interval k G [Mvr, (M + l)7r] and different M series do 
not overlap. 

In table ||, the bouncing ball eigenvalues (chosen by inspection of at the individual wave- 
functions) are compared with the quasi EBK quantization condition (l60| ) for a = 5. The EBK - 
solutions are also marked in Fig. ^. The importance of the additional terms fc/vr + 1/4 becomes 
evident when comparing the results with the spectrum of the rectangle obtained from to the 
condition a^/k'^/ir'^ - M'^ = N, i.e. k'^/ir'^ = + N^/a"^. A gain, we find that some bouncing 
ball states expected from the EBK formula ( |60D appear to have lost their regular nodal pattern 
due to interference with non-bouncing ball states. 

Note, that the T-operator for the rectangle is not recovered in the limit a ^ oo. The 



R matrix (56) is independent of a and the circle boundary can thus not be treated as a small 
perturbation even for large a values. The solutions of ( |60[ ) can approximately be written in the 
form 

^ = + ^{N - Ar,,M - ]f (61) 



with An,m = \/M2 + f j - M ^ for M » 



a 



A rectangular-like spectrum is achieved in the "integrable" limit o ^ oo, however, in a very 
special way. First of all, Eqn. ( ]6l| ) contains an extra phase 1/4 originating from the caustic 
in the circle, which is not present in the rectangle. In addition, the different M-series always 
have a finite cut-off for finite a, see (|60|). The number of states in a given M series increases to 
infinity only in the limit a ^ oo. 



The determinant det(l — T) is analytic in each strip Rek E]m7r, (m + l)vr[, with m integer, 
but has 8b cusp clt k = mvr. The non-analytic behavior is introduced through infinitely long tra- 
jectories contributing in leading order at kn ~ m, see Eqn. (]53|). These cusps are a consequence 
of omitting the bouncing ball contributions, which itself exhibit a cusp at integer multiples of 
vr, see Eqn. (^). The full spectral determinant D{k) ( |lO| ) is analytic and nonanalytic behavior 
in the bouncing ball contributions is canceled by the near-bouncing ball dynamics. 

Of special interest is, however, the maximal length of trajectories necessary to resolve the 
regular part of the quantum spectrum at a given wave number k. Inserting (^) in (^3|) leads to 
an estimate for the trajectory lengths contributing dominantly to the ground state (M, N = 1) 
in each M series, i.e. to the state located next to the cusp k = Mvr. We obtain 

4 2 4 2^ 

Lmax ~ rimax ^ -a M ^ -a -. (62) 
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Table 1: Eigenvalues kqm belonging to bouncing ball eigenstates compared with the EBK-like 
quantization condition (|60| ) for a = 5; here, A^, M denote the approximate quantum numbers; 
the numbers in brackets correspond to an enumeration of all states in successive order. 



This is in contrast to general semiclassical arguments for bound systems leading to a cut-off for 



periodic orbits sums at half the Heisenberg time |24|. This transforms for billiard into a cut-off 
in the length spectrum according to 

- A 

Lmax = ■ndik) ~ jA;, (63) 



where d{k) is the mean level density and A the area of the billiard. The estimate ( |62|) scales 
differently with the billiard parameter a and deviates thus especially in the "integrable" limit 



a — > oo. The cut-off (63) as derived in Ref. [24| is based on two main assumptions: the 
overall validity of the Gutzwiller periodic orbit formula and the analyticity of the semiclassical 
spectral determinant in a strip containing the real energy axis. Both assumptions which might 
be intimately related for bound systems fail here. Note, that this is not necessarily true for 
scattering systems. A semiclassical quantization of two examples showing intermittency, the 



Helium atom [35| and Hydrogen in a constant magnetic field for positive energy jSg], could 
be achieved within the Gutzwiller approach. However, the intermittent part of the dynamics 
introduces in these cases as well regular structures in the resonance spectrum and non-analyticity 
in the spectral determinant. 



4 The bouncing ball states 

The threshold (^) can be interpreted as the boundary of a region in phase space surrounding 
the marginal stable bouncing ball family. Semiclassical contributions from the dynamics in 
this region have a form similar to the one obtained for stable islands [||, ^ or in integrable 
systems The threshold values can be directly translated into momenta in phase space, i.e. 

p\. = ib/ssin^f with (pt ~ arctan(no,i/a) ~ arctan -^/fc/vr, see Eqn. (|l^). A semiclassically stable 
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island can thus be defined covering a phase space volume Vreg = k'^Vren with 



Vreg{k) = Aa d(/> 4a - - arctan W - (64) 




k/iT » 1. 



The size of the semiclassical island approaches zero in the limit k ^ oo (compared to the volume 
of the full phase space). The number of quantum states associated with this island thus increases 
like 



Nreg{k) ^ a +0{^k/7T). (65) 

on average. This estimate coincides with the average increase of regular states given by the 



EBK-quantization condition (^). It exceeds previous results by O'Connor and Heller |46|, 
who obtained a linear increase in k for the number of localized bouncing ball states up to the 
semiclassical limit A; — > oo. Note, that a fc^^^ increase in the number of bouncing ball states is 
still consistent with the Schnirelman theorem [^] , as the fraction of regular states compared to 
all eigenstates approaches zero in the semiclassical limit. 
The oscillating part of the level staircase function 

Noscik) = N{k) -N{k) (66) 

shows a strong periodic modulation in the stadium, see Fig. {N{k) = ^n^i^ ~ ^n) denotes 



here the quantum level staircase function and N its mean part given by the Weyl formula ) . 
This modulation coincides with the oscillating part of the bouncing ball contributions |21, ^% 



P3| |, see Eqn. (^. It was shown in section that the bouncing ball family does not contribute 
to individual eigenvalues. From the point of view of individual states in the spectrum, the 
oscillatory behavior in the level staircase function is caused by a periodic change in the spacings 
between neighboring eigenvalues. The spacings are unaffected by contributions coming from the 
classical bouncing ball family. 

Taking the concept of a semiclassical stable island seriously, we expect the spectrum of the 
stadium billiard to be divided into two different subspectra. The majority of states belong to 
the "chaotic" subspace formed by non-localized eigenstates (leaving aside the phenomenon of 
scaring along short unstable trajectories). Their number increases on average like 

_ _ _ A /fc\3/2 

N,haos{k)=N{k)-Nreg{k)^—k^-ai^-j + 0{k), (67) 

and the volume of the billiard is A = a + 7r/4. The level statistics of this subspectrum is 
expected to follow the GOE prediction of random matrix theory. The statistics as well as the 
average level repulsion is then stationary, i.e. independent of k, in the unfolded spectrum. The 
level staircase function is expected to be structureless showing only "statistical" fluctuations 
around the mean value Nchaos{k)- The regular subspectrum contains the eigenstates originating 
from a quantization of the bouncing ball island (|6^. The coupling of the dynamics in the 
"semiclassically integrable" region to the outer classical motion is weak compared to the mixing 
in the outer region itself. The level repulsion between bouncing ball states and chaotic eigenstates 
is thus small compared to the coupling among non-bouncing ball states itself. A possible 
structure in the level staircase function for the regular states which can be obtained from the 
EBK- quantization condition ( |60| ) can therefore survive in the full spectrum. This is indeed the 
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Figure 6: The oscillating part of the level staircase function (a) for the full spectrum; (b) for 
the bouncing ball states only. The full line corresponds to the oscillating part of the bouncing 
ball contributions, see @). 



level staircase function can be defined with the help of 



(68) 



case. The oscillating part of the EBK 
Eqn. (H), i.e. 

K!^{k) = NEBK{k) - Nregik). 

Here, NEBK{k) is given by the number of levels obtained from the quantization condition ( |60|) 
up to a certain k - value. The result is shown in Fig. ^3, (where the next to leading terms in 
Nreg havc been fitted numerically). The function N^g^ shows exactly the same modulations as 
the full spectrum and coincides also with the oscillating part of the bouncing ball contributions. 
We conclude, that the dominant oscillation in the level staircase function of the full spectrum is 
caused by a modulation in the density of bouncing ball states only. 



5 Conclusion 

I have shown that the marginally stable bouncing ball family in the stadium billiard does not 
contribute to individual quantum eigenvalues. The result is confirmed by a semiclassical quan- 
tization of the quarter stadium using Bogomolny's transfer matrix technique in a representation 
which excludes the bouncing ball orbits explicitly. The regular bouncing ball quantum states 
in the spectrum derive semiclassically from the near bouncing ball dynamics alone, and these 
states follow a simple quantization rule. The trace of the transfer operator can be approximated 
by periodic orbit contributions, which, however, show a Ti-dependent transition from Gutzwiller 
to Berry-Tab or-like behavior when approaching the bouncing ball family. This leads to the 
concept of a semiclassical island of stability surrounding the marginal stable family in phase 
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space. The boundary of this region is exphcitly ^-dependent and the phase space volume of 
the island shrinks to zero (compared to the total volume) in the semiclassical limes ?i — > (or 
k oo). The quasi EBK - quantization formula can be associated with a quantization of 
the semiclassical stable island. The periodic modulations in the level staircase function can be 
related to a periodic change in the density of bouncing ball eigenstates. 

The results demonstrate that averaged dynamical properties like ergodicity and positive 
Liapunov exponent are not sufficient to ensure the applicability of the Gutzwiller trace formula. 
The "chaoticity" of the stadium is in a semiclassical sense indeed ^-dependent, and the billiard 
is for small k closer to an integrable system than to a hard chaos one. 

The Gutzwiller periodic orbit weights have to be modified in the whispering gallery limit as 
well, see appendix An accumulation of periodic orbits towards a limiting cycle of finite length 
has also been found in other systems as e.g. in the cardioid billiard ||4^, |4^], in the wedge billiard 
[ pO| , and in the Anisotropic Kepler Problem |5^. We expect that a careful treatment as 
outlined in the appendix will solve problems concerning a semiclassical quantization of these 
systems. 



The results obtained here for the stadium billiard are expected to be generic for systems with 
mixed phase space structure (however complicated by the existence of island chains and can-tori 
surrounding the stable island itself). The stable island influences the classical dynamics in the 
outer chaotic region by creating intermittency. The regular regions appear to be larger than the 
actual size of the stable island due to the finite phase space resolution of quantum mechanics. 
The findings explain in a natural way the existence of EBK quantum states associated with 
regions outside a stable island as found in |5|, Q . The width of this semiclassically integrable 
layer is h dependent. A semiclassical quantization of stable islands as well as the behavior 
of localized wave functions on classical boundaries and the description of tunneling through 
dynamical separatrices ||5|, 12, ^ will be sensitive to this behavior. The results derived here 
indicate a failure of the Berry-Keating periodic orbit resummation [24| for mixed systems due 
to the intermittency introduced by the stable regions. 

The work presented here is restricted to semiclassical aspects. The influence of intermit- 
tency on Frobenius-Perron and related classical operators |51, 52] is so far best described by 
a so called BER-approximation |53|. The spectra of classical operators are directly related to 
classical [54, ^] and semiclassical sum rules as well as to spectral statistics |5^, 58| in 



hard chaos systems. The influence of intermittency on these results is still an open issue. 



Remark: 

When completing this article, I became aware of a recent work carried out by Primack et. al. 
[p9|]. The authors discus diffraction in the Sinai billiard by analyzing the Fourier transforma- 
tion of the true quantum spectrum in detail. They could indeed relate all deviations from the 
Gutzwiller trace formula to diffraction effects (due to the concave boundaries in this billiard) 
except for some of the near bouncing ball orbits. The results in |59] indicate clearly, that the 
influence of intermittency as dicussed in this article can also be seen directly in the Fourier 
spectrum. 
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A Monodromy matrix in the stadium biUiard 




The Monodromy matrix used in section y describes the hnearized motion near a classical trajec- 
tory in phase space in a plane perpendicular to the orbit and on the energy manifold. For billiard 
systems, this plane is spanned by the local displacement vectors 6p± pointing perpendicular 
to the actual momentum of the trajectory. An initial displacement is thus propagated according 
to 

where M depends on the path of the underlying trajectory. For 2-dimensional billiards, one 
obtains 

where t = L/k denotes the time between two bounces at the boundary, L is the length of the 
path and k = The contribution from a reflection at the boundary is 

M. = 




where k denotes the local curvature of the boundary, and Q is the angle between the orbit and 
the normal to the boundary. The k - dependence is scale invariant due to the transformation 
bqijz) = 5q{k = 1), 5p{h) = k5p{k = 1) and we may set A; = 1 in what follows. 

Bounces on straight lines in the stadium billiard are treated as free flights. The angle of 
incidence 6 is the same for all successive bounces of an orbit in the circle. The length between 
two of these bounces is L = 2b cos 9 where b denotes the radius of the circle. The Monodromy 
matrix for / free flights in the circle interrupted by (/ + 1) successive reflections on the circle 
boundary is thus given by 

/ -(1 + 20 2lbcose \ , , 

B The whispering gallery limit 

The whispering gallery limit of the Poincare map (|l7| ) is formed by trajectories (m, /) with m 
fixed and / — > oo in the notation of section |3.2| . These orbits approach the boundary of the 
billiard with an increasing number of bounces in the circle. Only m = and m = 1 is possible 
in the limit / ^ cxd. I will consider here the case m = and 6=1. The approximate whispering 
gallery map for m = 1 follows by analogy. 

There is one fixed point for each symbol pair (m = 0, /) with starting conditions 

= cos f Jt^) ; = 0. (70) 



,2 / + 1 

The angle of incidence for reflections in the circle is 

= (71) 
1 2 ^ ^ 
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The length of the corresponding periodic orbit approaches a constant, i.e. 

Li = 2a + 2{l + 1) cos Oi = 2a + tt - ^jf^ + 0{l~^). (72) 

The Monodromy matrix along the orbits starting on the Poincare surface of section can be 
deduced from appendix ^ and is 

( I 2a \ 

= I 2(;+l) 4a(;+l) (73) 
\ COS Oi COS Oi ) 

The periodic orbit weighting factor | det(l — Mi)| is thus 

I det(l - MOI = 4a-^ = -{I + 1)^ + + Oil-''). (74) 

cos 01 TT 6 

The determinant and thus the largest eigenvalue A; of the Monodromy matrix increase quadrat- 
ically with the symbol index. The main difference compared to the bouncing ball limit is the 
convergence of the period Li towards a finite value: the length of the billiard boundary. The 
whispering gallery limit therefore introduces no intermittency. The Liapunov - exponent of 
periodic orbits, 

A; = — ~ log /, 

M 

diverges logarithmically and the boundary orbit (0, 1 = oo) is infinitely unstable. 

I will show below that the isolated orbit condition is again violated when taking the trace of 
the transfer operator. A naive summation of Gutzwiller periodic orbit weights would overesti- 
mate the whispering gallery contributions considerably. The breakdown of the stationary phase 
condition is caused here by the accumulation of periodic orbits, i.e. stationary phase points, near 
the billiard boundary. 

The phase space area of starting conditions on the Poincare surface for trajectories with the 
same /-value is centered on the corresponding periodic orbit. Its size shrinks to zero in the limit 
I — > oo both in the x and 4> coordinate. The transfer operator in the whispering gallery limit is 
thus directly given by the Jacobian matrix (73) in local coordinates 6x = x — xi, 5(j) = (j) — (pi 
centered on the fixed points (70). The orbit length is approximately given by the quadratic form 

Li{5x, 5x ) = 

with 




T, — I dxdx dxdx' \ — I 2a 

' ~ d^L d^L - _J_ J_ 

\ dx'dx dx'dx' / x=x'=xi \ 2a 2a cos0; 




where L; denotes the length of the periodic orbit. When taking the trace of T we are interested 
in trajectories satisfying 6x = 5x' which demands (5(/> = 0. More important is, however, the 
restrictions on the Sx - interval, i.e. 

_ cos 01 cos 01 

2 (Z + l)(2/ + l) - 2 (2Z + 3)(2/ + l)' ^ ' 

The interval length and thus the integration range decreases like 
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TVT^^ = V -l=J^e^'^L,-qii+i)n 2 /'dxe*'^"', (76) 



The whispering gallery contribution to the trace of the T - operator are 
/27ri V 2a h 



where the upper and lower limits in ( |75| ) have been approximated by their mean value 
^ cos 6*/ vr^ 1 



TT 



(2/ + 1)(2/ + 3) ^ 2 (/ + 1)(2/ + 1)(2/ + 3) ' 



The interval length decreases faster than (the square root of) the prefactor in the exponent for 
increasing /. A stationary phase approximation is not justified for I — > oo, and we obtain (as in 
the bouncing ball limit) a transition behavior 

Note, that the number of periodic orbits which can be treated as isolated stationary phase points, 
increases only at a rate proportional to The contributions for large / fall off like and 

thus faster than in the bouncing ball case. In addition, there is a cancelation between the (0, V) 
and the (1, X) family which show the same overall behavior but appear with opposite signs. The 
whispering gallery contributions are thus small compared to near bouncing ball contributions, 
at least for a > 1. 
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